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Abstract: The "dialogue of multipoles" matched asymptotic expansion for small black 
holes in the presence of compact dimensions is extended to the Post-Newtonian order for 
arbitrary dimensions. Divergences are identified and are regularized through the matching 
constants, a method valid to all orders and known as Hadamard's partie finie. It is closely 
related to "subtraction of self-interaction" and shows similarities with the regularization of 
quantum field theories. The black hole's mass and tension (and the "black hole Archimedes 
effect") are obtained explicitly at this order, and a Newtonian derivation for the leading 
term in the tension is demonstrated. Implications for the phase diagram are analyzed, 
finding agreement with numerical results and extrapolation shows hints for Sorkin's critical 
dimension - a dimension where the transition turns second order. 
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1. Introduction 

The black-string black-hole transition (which includes the Gregory-Laflamme instability 
Q) is a phase transition in General Relativity which occurs at higher dimensions d > 4 
where black hole uniqueness fails. It is known to raise deep issues including topology change 
and critical dimensions Q as well as naked singularities and thunderbolts (see the reviews 
H [| and references therein). 

In order to test these issues, black object solutions were sought both analytically and 
numerically in l^xS 1 . One limit which is amenable to analytic study is the small black 
hole limit, where po, the black hole radius, is much smaller than L, the size of the compact 
dimension. Admittedly, small black holes are far from being the large black holes (po ~ L) 
which are involved in the phase transition and the latter should be studied numerically. 
Yet, they are still very useful in providing the following 

• Hints about large black holes via extrapolation. 

• Important tests for numerics. 
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In H we introduced a general perturbation theory for small black holes, an imple- 
mentation of the method of "matched asymptotic expansion" to the static case, which was 
termed "a dialogue of multipoles" since it describes how field multipoles change the black 
hole's shape (or mass multipoles) and these in turn change the field and so on. This method 
requires to obtain a solution in two zones: the near horizon zone and the asymptotic zone, 
solutions which must agree over the overlap region. 

In principle it should be possible to devise a perturbation theory in a single zone. For 
instance, Harmark wrote an approximate solution in a single zone ||, but nobody extended 
it to a full perturbation series. However, it is probably simpler to use the two-zone method 
since it does not involve the arbitrariness of an initial guess, and since the differential 
operator to be inverted is the same at each order. 

In H we obtained the leading order corrections to the metric. Recently a "matched 
asymptotic expansion" was used to obtain the explicit solutions in 5d up to the next to 
leading order ( O(p 4 ), 0(L~ 4 ) ) |7]] (see also || for a closely related work in the braneworld 
context). 

In this paper we extend our method to the next to leading order in the asymptotic zone, 
namely to Post-Newtonian order, for arbitrary d (namely, order 2(d — 3)). It turns out that 
at this order a new qualitative phenomenon appears - divergences. Such divergences are 
familiar from Post-Newtonian studies, but we are not aware of a general prescription that 
works at all orders. In section || we obtain such a clear regularization in terms of a "cut- 
off and match" method, which relies on Hadamard's "Partie finie" regularization and is 
summarized in equation ( |2.15| ), This method shows some similarity with the regularization 
of quantum field theories (see Q for related ideas). Alternatively, a certain subtraction 
of self-interaction is shown to be an equivalent method at this order. We believe this 
equivalence should continue to hold at arbitrary orders, but this is not manifest in the 
current setting. Moreover, we find that at this order the matching is quite trivial and 
matching constants vanish. Again we do not know whether this can be made to persist to 
higher order and whether it depends on a clever choice of gauge order by order. 

The above-mentioned divergences did not show up in neither Q nor JJJ: || was using 
a different, single patch, method while in J?]] the authors succeeded in obtaining full explicit 
solutions for the metric without the use of Green's function, thereby circumventing the issue 
of divergences. Still, the regularization of these divergences is important both conceptually 
and practically: conceptually, they are a general feature of the two-zone method and 
show up at all high enough orders; and practically, they are essential if one wishes to 
use Green's function either while resorting to a numerical solutions for the whole metric, 
or when computing analytically the solutions' asymptotics, namely the thermodynamics. 
Moreover, similar divergences are very common in Post-Newtonian studies in general. 

The pre-existing results for the leading corrections are reviewed in section ^. Proceed- 
ing to next to leading order, we met difficulties in trying to obtain explicit expressions for 
the full metric, but we were able to compute the next to leading corrections to the black 
hole's mass and tension, using our regularization procedure, thereby confirming the results 
of @ with our different method. The computation is described in section || and the results 
are presented in subsection The leading behavior of the tension is given an intuitive 
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Newtonian explanation in subsection 4.4 . A higher order expression for the "black hole 
Archimedes effect" is given in Appendix [B|. 

Finally, we discuss the implications for the phase diagram in section [|, in light of 
the objectives presented in the beginning of this introduction. Within the small black hole 



region of validity we find very good agreement with numerical data in 5d and 6d jl(J, 1 1 , 0] . 
Extrapolating beyond this region we find evidence for a second order phase transition at 
large enough dimensions, which is an indication for Sorkin's critical dimension |IT 



1.1 Set-up 

Figure [l] illustrates the coordinates which we use. We denote by z, r the "cylindrical" 
coordinates, where z is the coordinate along the compact dimension whose period we 
denote by L, and r is the radial coordinate in the extended M. d ~ 2 spatial dimensions. In 
addition we introduce "spherical" coordinates p, x- Po is a characteristic size of the horizon. 
We are interested in obtaining solutions for small black holes, po <C L, so po/L is the small 
parameter. The "spherical" coordinates will be more natural in the vicinity of the black 
hole since it is nearly spherical in the small black- hole limit. 

A 



z 




- r 



V 

Figure 1: Illustration of the (r, z) "cylindrical" coordinates and the (p,x) "spherical" near horizon 
coordinates. The period of the compact dimension (in the z direction) is denoted by L and the 
radius of the black hole is denoted by po- 



The general method of "dialogue of multipoles" matched asymptotic expansion was 
explained in ||. The static nature of the problem under study makes this application of 
matched asymptotic expansion more transparent than the usual 4d applications. In this 
method we consider two zones (see figure ||): 

• The asymptotic zone where p 3> po and po is the small parameter. The zeroth order 
solution is flat space with a periodic coordinate z ~ z + L and the point at the origin 
(r, z) = (0, 0) removed. 

• The near zone where p <C L and L^ 1 is the small parameter. The zeroth order 
solution is the Schwarzschild black hole with radius po- 

Thus in the asymptotic zone we expand the metric in po and in the near zone we expand 
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I i y 1 1 re 2: The division of the spacetime into two zones: the near zone p -C L where po is 
fixed and the perturbative parameter is L -1 , and the asymptotic zone p ^> po where L is fixed 
and the perturbative parameter is pq. The two zones overlap over the overlap region which 
increases indefinitely in the small black hole limit. During the perturbation process the two zones 
are separate, and communicate only through the matching "dialogue" . The near zone is defined by 
{(l°> X) '■ P> Po} while the asymptotic zone is defined by {(r, z) : r > 0, z ~ z + L}\(0, 0). 

it in L- 1 

oo 

^ Symp) (x) = ^ Po ^r ymPj) (-) 5 

j=0 

oo 

3=0 

Orders which are integral multiples of d — 3 have a special role, and therefore we introduce 
a special square brackets notation for them 

9tl-=9%^. (1-2) 

The metric in the two regions must be consistent over the overlap region po <C p <C 
L, in a double expansion in po, L. This is the "matching procedure", whose first steps 
are summarized in figure ||. For example, the first match is between the Schwarzschild- 
Tangherlini solution in the near zone (order L~°) and the Newtonian approximation (order 

2. Divergences and regularization 

Divergences. When we proceed in the "dialogue of multipoles" matched asymptotic 
expansion beyond the leading order we encounter a new qualitative feature - divergences. 
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These occur in both the asymptotic zone and the near zone, and for concreteness we shall 
first describe the divergences and their regularization for the asymptotic zone in detail, 
and later the near zone will be described briefly. 

At each order, k, in the perturbation series we get an equation of the form 



L h (k) = Src( k) 



(2.1) 



where L is the linear operator which appears at first order, is a function which deter- 
mines the perturbation to the metric at order k and Src^ is a source term which depends 
on metric functions from lower orders. More specifically, in the asymptotic zone we have 



A h {k) = Src {k) ~ 8 kl $> k2 , 
where A is the Laplacian, $ is the Newtonian potential given explicitly by 

1 



d-3 ' ' 

(r 2 + (z + nL) 2 ) 2 



(2.2) 



(2.3) 



d kl is a symbolic notation for some differential operator of order k\ and k\,k% > 2. 1 
For small p behaves as 3> ~ l/p d ~ 3 and for large r as $ ~ l/r d ~ 4 . For example the 
Post-Newtonian equations (PN) are 



A h, 



m 

y 




, 



d ~ 2 * * 

<t> , <J> a 

2{d - 3) J ' 



(2.4) 



where i,j run over the spatial indices. 

Usually we solve an equation such as (2.2) by means of a Green's function. Namely, 
we take the solution to be 



hc(x) := / dx' G(x,x') Src(x') 



where in this case 2 



G(x, x) oc 3>(x — x ) 



(2.5) 



(2.6) 



Divergences of the integral may come from regions where the integrand is infinite or 
from large regions of integration. In our case short-distance divergences are possible in the 
vicinity of x' — > (there are none as x' — ► x), and long-distance divergences are possible as 
r — > oo. By a standard abuse of language we shall refer to the short-distance divergences as 



1 This form allows for non-linearities arising from both the Newtonian approximation ~ $ and sources 
introduced through multipole matching ~ d l &. k,2 > 2 since the source is non-linear and fci > 2 from a 
dimensional analysis. 

2 We take the Green function with no boundary conditions at 0. One could have taken a more involved 
G(x, x') with an arbitrary multipole distribution at x — 0, but this would not cure the divergences which 
we will encounter. For the precise prefactor in (2.6) see (4.14). 
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"UV divergences" and to long-distance ones as IR even though the theory is not quantum 
and distances and energies are not reciprocal. 

UV divergences. The integrand of fl2.5|) diverges as p' — > due to the singularity 
in Src(p') which is inherited from <£(//)• Assuming that for small p' the behavior of the 
singular and angle-independent ("S-wave") part of the integrand is G(x, x') Src(x') ~ p'~ s 
we get from ( ^5| ) 

ha(x)~ I dp 1 P ' d - 2 P '- s . (2.7) 
Jo 

Therefore we have UV divergences (for all x) exactly if s > d — 3. For PN we see from 
Q2.4D that Src(p') ~ 3?^ &j ~ p~ 2 ( d ~ 2 \ 3 and therefore we indeed have UV divergences for 
all d. At higher orders the divergences only get worse. 

IR divergences. If for large r Src behaves as Src(r) ~ r~ s °° we get 

/OO -I 
dr'r' d ~ 3 r'~ Soa . (2.8) 

Therefore we would have divergences exactly if Sqo < 2, but the slowest relevant decay rate 
of Src - a PN source in d = 5 - is Sqq = 2{d — 2>) = 4 and hence there are no IR divergences. 



Unregulated multipoles. Even though the Green function integral diverges we can still 



find a particular solution of (2.2). Locally around the origin we may separate the equation 



by using radial variables (p, %) and we find 

g; + ^l df Jl±lzR) hl = Srci . (2 . 9) 

p p p l ) 

The particular solution is 

h oc p- Sl+2 , (2.10) 

where the constants s\ are defined through Srci(p) ~ p~ Sl , and one should allow for a 
possible addition of a homogeneous solution of the form p l , p~( d ~ 3 )~ l . 

Indeed, if particular solutions can be found in explicit form, and their homogeneous 
part determined by matching, then there is no problem of divergences, as was done in [f?j for 
d = 5. However, we are often interested only in the asymptotic form of the metric and then 
a Green-function-based multipole expansion is useful. For instance, in this paper we will 
compute only the corrections to the thermodynamic quantities, the mass and the tension, 
rather than the whole metric. Rephrasing from a different perspective, the particular 



solution (2.10) at small p does not give information about the large p asymptotics, and the 
separation of variables is not valid globally. 

We now wish to re- formulate the divergence problem in terms of multipoles, which are 
particularly suited for divergences localized at a point. Since the divergences come from 
the vicinity of x' = we Taylor expand G(x, x') there. The convergence radius p cnV g would 



3 More precisely, one needs to find the degree of the singular and angle-independent part of 
G(x,x') Src(x'). Since for PN the leading singularity of Src(p') has a quadrupole part then s > 2(d — 3) 
and the conclusion is unchanged. 
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be ~ min(L, x) and for x — * oo p cnvg ~ L. Then we separate the domain of integration 
in ( |2,5[ ) to within the convergence radius p < p cnvg and outside it. Outside there are no 
divergences, and inside we may perform a standard multipole expansion 

h{x)^d' h ...d' i[ G(x,0) f dV Src{x')x' h ...x' il . (2.11) 

J P^pcnvg 

Thus we see that all the divergences may be encoded by the diverging multipoles at the 
origin, defined by the small e behavior of 

Mr.= I dV Src(x) x h . . . x H . (2.12) 

J\p\>* 



Regularization. The matching boundary conditions (b.c.) suggest a natural regulariza- 
tion procedure: "cut-off and match". Namely, cut-off the integral ( |2.12| ) at p = e, adjust 
it by a constant (corresponding to adding an allowed homogeneous piece to h) to conform 
with the matching b.c, and finally take the limit e — > 0. 

This idea is translated into formulae as follows. Expand the multipoles in a Laurent 
series in e 



M l {e)= ^ Mje J , (2-13) 

3 = 3min 

where at order k dimensional analysis gives j m i n = — (k — I — (d — 3)). From the finite 
part multipoles, o, we may construct a function h which is known in the mathematical 
literature as Hadamard's partie finie [14], or Hadamard's finite part, and hence we shall 



denote 

Pf(Mi) := Mi,q (2.14) 

(see also |l5|] for its use in the Post-Newtonian context). Alternatively, the regularization 
could be defined through analytic continuation and it is also closely related to Cauchy's 
Principal Value. 

The divergent piece in ( 2,13] ), (namely the sum over j < 0), is interpreted as coming 



from the vicinity of the origin, and thus according to the matching b.c. should be replaced 
by an appropriate value of the multipole - A4J na * c/l read from the asymptotics (p 3> po) of 
the near zone metric. Thus we arrive at the expression for the regularized multipoles which 
is one of our central formulae 



Mi -»■ M t R = Pf(M) + M x atch 



(2.15) 



"Renormalization" . This regularization procedure is reminiscent of renormalization in 
Quantum Field Theory. There we know that in order to renormalize a theory it must have 
as many free parameters, or counter-terms, to be set by experimental data, as divergences. 

4 In QFT a theory is called normalizable if the number of counter-terms is finite, and non-renormalizable 
otherwise. In our case there will be infinitely many divergences, but all the necessary "experimental data" 
will come from matching. 
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Here by analogy we need 

# (divergences) = #(matchings) . (2-16) 

However, at each order there are only certain (finitely many) matchings coming from 
the other zone, as encoded by the "perturbation ladder" (see figure ||). What would 
guarantee that ( 2.16| ) is satisfied? The answer turns out to be reminiscent of Quantum 



Field Theory as well - it is dimensional analysis! At order k in the asymptotic zone we are 
considering divergences of the form 

Eo. (o 17) 

p i+d-z e a L b ' \*- Ll ) 

where a > 0, b > 0. Since all the quantities in this expression have dimensions of length, 
but the expression itself (namely h) is dimensionless then a + b = k — I — (d — 3) and in 
particular k— id— 3) — I > 0. However, this is exactly the condition for matching, as is seen 
by considering appropriate terms at order b in the near zone, namely taking a = 0, b > 0. 

No self-interaction. We are also familiar with a different approach to regularization 
- the removal of self-interaction. The idea is to separate the black hole solution into a 
Schwarzschild piece and a correction coming from the compact dimensions. Intuitively one 
expects that all the divergences come from the singular Schwarzschild part, and that these 
will be absorbed automatically by adjusting its parameters. 

This approach is commonly used in the Post-Newtonian context (see for example the 
review talk |16[| ). However, it is not known how to implement this idea at arbitrary orders 
of perturbation, but rather some results are known at low orders. We expect the "no-self- 
interaction" to be equivalent to our method and it would be interesting to demonstrate it. 
Our explicit calculations in the next chapters do lend support to this idea, as it turns out 
that if one omits singular products from the double sum over images in PN, it is equivalent 
to Hadamard's regularization. In other words, Hadamard's partie finie of the singular part 
is zero at this order. 

Moreover, we observe another simplification whereby the matching is exactly cancelled 



by the multipole contribution of the 0(<I> 2 ) term within the Laplacian in (2.4). This may 
very well be a special feature of our gauge, and it would be interesting to determine whether 
it can be made to continue in higher orders as well. So altogether in our PN computation 
regularization may be replaced by discarding the singular self-interaction source term and 
matching is automatically cancelled. 

The near zone. In the near zone the regularization of divergences works in parallel to 
the method just described for the asymptotic zone, so we only describe the main points. 
Actually, in this paper we shall not treat the near zone at all. 



Here the linear operator, L, in (2.1) is given by the Heun equation (see |5[), which 
for /) > po reduces once more to the Laplacian. The divergences are IR rather than 
UV. Regularization proceeds by introducing a large-distance cut-off R and considering the 
Laurent series for multipoles as R — > do. After a partie finie regularization the matched 
value from the asymptotic zone should be added just as in ( p. 15 ). 
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3. Review of leading order results 



3.1 The general procedure 

Using the matching procedure described in || we obtained the leading corrections to the 
Schwarzschild metric of a black hole as a result of the compact dimension. The matching 
procedure for the leading corrections is summarized in figure ||. 



Newtonian 
Approximation 
d-3 



PN 

Approximation 
2 (d - 3) 



Asymptotic Zone 




Figure 3: The first steps in the matching procedure. Each box in the top row denotes a specific 
order in the asymptotic zone, and the lower row depicts the near zone. Arrows denote the flow of 
matching information between the zones: Dark (red) arrows denote monopole matching (7 = 0), 
light arrows (light-blue) denote the quadrupole (1 — 2), and higher I arrows are not shown. 



The metric in the asymptotic zone is expanded in the radius of the black hole po in 
the following form 



(asymp) 



+ + h®, + 



(3-1) 



where h^l oc pQ~ 3 ,hfH oc p^ d 3 ^ and so on. Note that the matching procedure may 
introduce corrections which are not integer powers of p^~ 3 but this would happen only 
after the post-newtonian corrections. 

According to the matching procedure, in order to determine the Post-Newtonian cor- 
rections (order 2 (d — 3) in figure |||) we first need to know the Newtonian approximation 
and the near zone leading order. Therefore we start here by reviewing the results that were 
obtained in ||] for these orders. The leading order in the asymptotic zone is simply the 
Newtonian potential while the leading order corrections in the near zone are decomposed 
into multipolar static perturbations to the Schwarzschild- Tangherlini metric (following 



Regge- Wheeler [17] and using a similar gauge). 

For the post-Newtonian expansion we make a quite-standard gauge choice. One starts 



by writing the Ricci tensor as [18] 



R 



a/3 



d 2 g' 



pa 



■ 2 9^9u P g dxadxf3 



r 



(3.2) 



- 9 - 



where r„ a g and Tff are the Christoffel symbols of the first and the second kind, respec- 
tively, and in addition one defines 



v - l ( dgPa ^ 

^ - 2 I 9w> dx „ + 9ua dxfM 9w9vo dxa 



Next one chooses the Harmonic (or de Donder 5 ) gauge by the requirement that 

r" = Or" = J=JL(^g-gP><) = 0, (3.3) 

where we denote by g the determinant of the metric g^ v . In this gauge, the last term in 
the expression of the Ricci tensor above vanishes. 

3.2 Newtonian approximation 

The first step in the iterative procedure was to solve the equations in the Newtonian 
approximation (see for example |5|, ||, p0[ ) 

t& ] = *■ (3-4) 

and 

*--=pt 3 E — 71EF' ( 3 - 6 ) 

„=_oo (r 2 + {z + nL) 2 ) 2 

where the Latin indices stand for the spatial components. In order to match this ap- 
proximation with the near zone we have to expand around p = in polar coordinates 

(p,x) 



PA~* , o . ti a ^A d ~ 3 ) ( d ~ 2 ) P 2 pV 
P d 



^ z ) = ^3+ 2 -^T3C(d-3)+ ~> ' ((d - 1) U 2 /( x ) + 0{p% (3.7) 

where £ is Riemann's zeta function 6 and 



nf(x) (cos 2 ( X )(d-l)-l), 

which is the generalized Legendre function that corresponds to the quadrupole. In partic- 
ular, the overall prefactor in ( |3.6[ ) is set by the matching described by the left most arrow 
in figure || 

From the asymptotic form of the metric (monopole terms), when r 3> z, 

9^ = ^1 + (^3 ) , (3.8) 



3 The first introduction of this gauge appeared in Jis|. 
"Riemann's zeta function is defined as £(s) = X^^Li "tt- 
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one can extract the mass and tension ||i"o| , 21 denned through the standard first law 



dM = — dA + rdL . (3.< 



57T 



(For weak sources this definition is equivalent to M := J dV d -iT tt , r := — J dVd-iT zz /L.) 

In harmonic gauge the asymptotic measurable quantities can be expressed in terms of 
the constants in the asymptotics as follows [1C, 21 1 

M = -^r((d-3)c tt -c zz ), (3.10) 

(c tt - (d - 3)c zz ) , (3.11) 



16 vr G d L 
where 

2 d ~ 2 

Qd ~ 3 = i#) ' 

is the area of a unit S 3 . 

At this order, the asymptotic constants can be determined from the asymptotic ex- 
pansion 7 (r 3> z) of the Newtonian potential 

_ n d _ 2 d-3 1 . „„ / 1 

* ~ f 

Then 



* = s£w r ^ +0 l^<- (3 - 12) 



c [1] - 




2 




Vd- 


-3 


c [l] = 


Vd- 


-2 


^zz 




-3 



%r 3 , (3.i3) 



d-4 

^4"" <3 ' 14) 
Thus, the tension vanishes in the first order approximation and for the mass we obtain 

3.3 The near zone metric 

In the near zone we obtained the first monopole correction to the Schwarzschild metric by 
matching with the Newtonian approximation ||. The metric in the near zone reads 

n near _ ( pt 3 \ ( i 2gd-3) P d ~ 3 \ , x s 



,near 



i-o ii+ 2 -mi^)+o( T ^), (3.16) 



..near _ n 2 A , 2j(d-3) P d - 3 \ n / i\ 

The metric in the near zone is written in Schwarzschild coordinates when the metric 



in the asymptotic zone is written in Harmonic coordinates (ones which satisfy (3.3)). We 



Which may be gotten from flux conservation for 
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denote by p s the radial coordinate in Schwarzschild coordinates and by p the radial coor- 
dinate in the Harmonic coordinates. In order to match them we have to use the leading 
terms of the transformation rule || 

P = Ps- 9( /\, d - 4 +0(^- 3 )- (3.17) 
2 [d - 3) p s p s 

and get the near zone metric in Harmonic coordinates 

^ = j 1 _,g)( 1 _m^\ +0 (ii), ( , 18) 




p2(d-3) ' L d-2 

near _ near | , /'m 1 I , 2£(rf - 3) Pq 3 \ ^ 1 I 



(d-3)L d - S J \p2(d-3) ' L d ~ 2 

Note that the prefactor of in the near zone determines the monopole corrections to 
the asymptotic zone metric through the matching procedure. In particular, the term that 

2(d-3) 

behaves as ^ d ° 3 Ld _ 3 is the one that should be matched with the post-Newtonian corrections 
to the metric in the asymptotic zone (see figure |||) . 

Higher multipoles. 

From the "perturbation ladder" (figure ^) we see that proceeding in the near zone 
beyond the leading correction (order d — 3) we have contributions of higher multipoles 
that are determined from the Newtonian potential, before the first non-linear iteration 
at order 2(d — 3). The multipole moments of the newtonian potential are matched with 
the multipole linear perturbations of the BH. In our case, only even multipole numbers 
contribute corrections to the metric due to the x ~> n ~ X symmetry. Any even multipole 
number I < d — 3 contributes a correction to order I + d — 3 in L _1 before order 2 (d — 3) 
(see H), where the next monopole correction enters. For example, in d = 5 there are no 
corrections of higher multipoles before order 4, in d = 6 we have a quadrupole correction 
I = 2 and so on. 

For any / < d — 3 we obtained in || explicit expressions for the correction to the near 
zone metric. We review here the final results. The correction of order I + d — 3 can be 
expressed in the form 

5 (near,^-3) ^ ^ = g (W +d -3) ^ + ff (near,/ +d -3) dp 2 + ^(near^-3) ^2^ (3 1Q) 
where 

9xx ~ Ps c iPo -r^ + _2j_) _ E i( X ) U o (X). I 3 - 20 ) 

E A X) = (l - (d - 3) X A) ((1 - X) *<! - ^,2 + ^,2 ; 1 - X)) , 

d-3 

X := (3.21) 
Po 
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2-Fi is the hypergeometric function and H l Q d (x) are the generalized Legendre polynomials 
given by a Rodriguez formula 



sax (x) 



r(z + 



1) 



1 



d 



sin(x) dx 



sin 



2Z+d-4 



(x)- 



(3.22) 



The constants cj come from the metric in the asymptotic zone, namely, from the Newtonian 
potential. These are the constants that appear in 



(asymp,d-3) 
9xx 



d — 3 d 



(3.23) 



- see ( |3.7D for the explicit values of c% for I = 0, 2. The other two components of the metric 

,. , , • i r (near,z+d-3) 

correction are determined from g^ x 



9tt 



(near,/+rf-3) 



c l Po 



i r(i+^)r( 2 +^3) 



^f Al (p s )I^ d ( x ) 



r(i+^)(i-i) 



Z < d- 3 



(near,/+d-3) / g(i+jpsj r ( 2 +5^) B t (p a ) ttM 



(3.24) 



2pp 



r(i- 



0' 



where 



/ := 1 



d-3 ' 



(3.25) 



A(/0, 



(d-2)p s f 



s dEi d 2 Ei 



dp s 



21(1 + d- 3) 
and Bi is obtained from the algebraic relation 

B l = -A l -(d-A)E l . 



dpi 



+ 



d-4 



4. Post-Newtonian Thermodynamics 
4.1 Post-Newtonian equations 



Using the results for h\J, (the Newtonian approximation) we can write Einstein s equations 
for h^l, i.e. the Post-Newtonian equations, in the following form 

A (ft!? + ^* 2 ) = 0. (4.1) 

^ tj 2 l] \d-3J J 2(d-3) ' J 

The origin contains a singular source which will be accounted for by matching rather than 
by introducing a point-like source. 
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Let us look at the "source" term (the RHS) for the spatial components. This term is 
proportional to 

2 2(d-3) (xi + mL5 iz ) (xj + riL5 jz ) 



ez + (2 + mi) 2 ) 2 (r 2 + (z + nL) 2 ) 2 

In this sum we see the effect of interaction between the black hole and its mirror images, 
including the interaction between the mirror images themselves. One term in this sum is 
different - the term that corresponds to the self-interaction of the black hole, namely, when 
m = n = 0. This term exists in the post-Newtonian equations regardless of the boundary 
conditions. It has a singular behavior near the origin being proportional to s^Tj • For 
i = j its integral over the volume diverges like for e —* 0. 

We shall follow the "no self-interaction" approach together with "matching regulariza- 
tion" (see the discussion in section |2|). Since the equations are linear we can separate the 
equations to the self-interaction part (SI) and the regular part (REG) which consists of the 
interaction with the mirror images. In order to keep the periodic boundary conditions in 
the compact dimension, we include in the SI part the self-interaction terms of the mirror 
images even though they are not singular in the domain \z\ < L/2. 

Accordingly we separate 

tfl = h$ 81 + h^ G . (4.2) 

The separated PN equations are 

^(h^ SI + ^h)=0, (4.3) 



where 



2 



3,2 _ 2(d-3) V- i 



and 



a(/4?^ g + ^w)=o, 

2(d-3y dij ^ REG ) = ~2(d-3) 



(4.5) 
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where 



.2 _ 2(d-3) y- 
^ REG — Po 2^ 



jkn {r 2 + (z + m L) 2 ) 2 (r 2 + (z + nL) 2 ) 



d-3 ' 
2 



(d -3)^ (d - 3) E 



3 1 REG 



(xj + mL 5 iz ) (xj + nL 5 jz ) 

^ n (r 2 + (z + mL) 2 p (r 2 + (z + n L) 2 )^ ' 

( m, n E Z ) . 

Finally, the boundary conditions determining the homogeneous solution are supplied through 
matching with the near zone metric ( |3.18| ), as usual. 

The self interaction part — SI. An explicit form of a particular solution for the singular 
part (4.3,4.4) can be found. For this purpose let us look at the equations in the case of a 



[2]0 



single source without any mirror images (n = 0), whose metric we denote by h^ v 

■ 0, 



/ r , 2(d-3) ' 

A h [2]0 + -^ 

I tt + 2p2(d-3) 



A h 



2]0 



1 



Po 



2(d-3) 



2 (d - 3) 2 p 2 ( d ~ 3 ) 13 J 
The solution for the equations above is 



(d -2)(d- 3) 2 (d-3) XjXj 
2 P ° p 2(d-i) 



(4.7) 
(4.8) 



2(d-3) 

Po 

2p2(d-3)' 



(4.9) 



, [2]0 _ 2 (d-3) 
n ij — Po 



'ln(p) / 4xj rr 



t- \ , Xj 11 Sjj j — p- 

~ °U ) + 12^" 96^ 



2(d-3)V (d ~ 3) 



Xi Xj (d— 3) 



d > 5 



Now, the solution to the full SI equations ( pL3; , |4.4| ) can be obtained by taking z 
z + n L in (WM) and summing over n£Z 



n tt 



h [ 3 SI 



h [ S°(r,z + nL), 

=— oo 
oo 

£ ^ ] °(r,* + nL). 



(4.10) 
(4.11) 



The regular part — REG. The particular solution in the asymptotic zone can be ex- 
pressed by an integral over Green's function 

f% ]REG = -\*%EG, (4-12) 

^ REG = 2-{^ 5 ^ REG ~ 2jd^j I (^)REG(r'^')G(r-r',z-z')dV d ^, 

(4.13) 
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ht 


hzz 


SI 








REG $ 2 shift 


— 2 C$ + 


2(d-3) 2 C * 


REG Green's 








order [1] near 


-2C 


, 2C 



Table 1: Matching budget for both h u and h zz at Post-Newtonian order, namely for the coefficients 
of O ^Po^ _3 V(/ 0<i_3 L d ~ 3 )j . At PN these quantities have 3 contributions SI - Self-interaction 
solutions see (4.9); REG <f> 2 shift - namely the contribution from &% EG to the regular term see 
(4.12,4.13,4.17). Here c$ := 4£ is the coefficient of the relevant term in Q 2 REG (4.17) and in this 
table C : = Q{d — 3); REG Green's - for h zz the Green's function integral in (4.13) has no relevant 
term according to (4.18). For htt this term is identically zero (4.12). In the fourth line we see the 
leading order results for the near zone, taken from (3.17), and we confirm that they are identical 
to the sum of the 3 items above them and therefore the matching constants vanish. 



where Green's function is 

oo 

G(r-r',z-z') = - £ ^. (4.14) 

(d - 3)Oi_2 n f^oo ((r - r') 2 + (z-z' + n L) 2 )~ 

4.2 Matching 

We still need to match the solutions with the near zone, a matching which will turn out to 
vanish. As encoded in the "ladder" figure (figure ||) we should match the monopole term 
at order 2(d — 3) in the asymptotic zone, with the monopole of order d — 3 in the near zone. 
Accordingly, the relevant term in the double expansion in the overlap region po <C p <C L 
is 

2(d-3) 

Po (4.15) 



pd-3 ]^d-3 ' 

The matching "budget" is composed of several items summarized in the table ||, which 
we proceed to explain. 

• SI — In the case of the Self-interaction equations we have an explicit solution (|4.S|). 



[21 [21 

One sees that neither h tt nor h z £ contain any term of the form ( |4.15| ). 



In the regular solutions ( 4.12| .4.13), we distinguish two parts, the first term is just a 



shift proportional to <&\eq, which we refer to as "REG < I )2 shift", and the second is 
an integral with a Green function kernel (there is no such term in ( 4.12|) ), which we 
refer to as "REG Green's" . 

REG <]? 2 shift — In ^\eg we nave a term of the relevant form ( [4 . 1 5| ) 
^REG = C(d " 3) + (t\d - 3) - \ C (2 - (d - 3))) + (4.17) 
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whose coefficient we denote by c$ = 4£(<i — 3), where we used the expansion for 
( |3.7| ). We write here explicitly only the two monopole terms, i.e. terms which are 
independent of the angle x- 



REG Green's — Next we show that the piece of h zz in (4.13) involving an integral 
with Green's function 



does not contain terms of the form ( 4.15| ). For this purpose let us look at the leading 
term near p = in the integral 

poo roc ii i 

I (<S>^, 3 ) REG ( P ')G(p,p')p' d -*dp'^ I (4.18) 

where in the first equality we estimated the singular angle- independent part of (&,i&,j)jiEG 
by l/p' d ~ 3 and replaced l/|x — x'\ d ~ 3 in the Green function by l/p' d ~ 3 with a cutoff 
at p. Thus, we find that the leading behavior when p — > is O (p~^ d ~ 5 ^), namely 
there are no terms that behave as ( f4,15| ). 

[2] [2] 

Adding up these 3 contributions for both h tt and h zz we find exact agreement with the 
near zone expressions ( |3.18| ), as seen in table [l|. Therefore the matching constants indeed 
vanish. 

4.3 Results 



We did not find an analytic solution for the integral ( 4.13 ), and hence we do not have a 
full analytic form for the metric at this order. However, we can analytically evaluate the 
asymptotic form of this metric, thus determining the PN corrections to the thermodynamic 
quantities, M, r. The results agree with Harmark's results || which were obtained using 
a different, single patch, method. 

[2| [2] 

Our goal is to find c tt and c zz , the next to leading corrections to gu and g zz for r>z. 
From the asymptotic form ( |3.12| ) we see that the terms in the expansion of $ 2 for r > 2 do 
not contribute to the asymptotic constants. Moreover, we found in the last subsection that 
the matching constants vanish. Thus only the source terms (the RHS in ( [4.1|) ) contribute 
and we have 

eg 1 = 0, (4.19) 

Pf stands for the finite part of the divergent integral according to Hadamard's partie time 
regularization. 

On dimensional grounds 



ft \ 2{d ~ 3) 

Pf f / {$ >z fdV d -i \ = const , (4.21) 
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where const = const{d). We evaluated these constants by separating the integral into 
SI and REG pieces. The partie finie of the SI part vanishes and we find (the detailed 
calculation of the integral and its regularization appears in Appendix [A]) 



pf / {^ )Z ydv d ^ 



2 (d-3) 

Po 

L d-3 



O d _ 2 (d-4)(d-3)C(d-3) 



Substituting back into ( 4.20 ) gives us the final expressions 



J2] 



0, 



Po 



2 (d - 3) C(d-3)(d-2)O d _ 



L d-3 



2Q« 



(4.22) 

(4.23) 
(4.24) 



Using ( |3.10 , |3.li|) and adding it to the results in the Newtonian approximation ( |3.15 ) we 
obtain the "measurable" quantities as 



M 
tL 



(rf-2)n d _ 2 d-3 

16tt G d PO 



1 + 



C(rf~3) £t 



L d-3 



+ o[p 3 {d - 3) 



i-3) (d-2)fl d _ 2 gd-3) pl {d 3) n ( 3(d-3) 
L d- 3 -\-U\p Q 



32tt G d 



(4.25) 



4.4 Newtonian derivation of tension to leading order 

We shall now show that the leading order tension ( 4.25| ) can be understood from Newtonian 
gravity. From the first law of black hole thermodynamics ( |3.9| ) we can express the tension 
as 

(dM\ 

t = ( — 1 , (4.26) 



dL 

namely, the change of the mass in reaction to a change in the period of the compact 
dimension when we keep the entropy constant. The demand for constant entropy of a 
small black hole just means that po (as well as the rest-mass Mo) are fixed. The black hole 
mass in the Newtonian limit is 

M = M + U , (4.27) 

where the mass is expressed as a sum of Mq, the intrinsic mass of the black hole given by 
the leading order of ( 4.250 

(d- 2)Q d ^ 2 d-3 (A 9S nj 

U <C Mq is the gravitational potential energy, written as a sum over all mirror images with 
the standard prefactor 1/2 to avoid over-counting 



U 



(4.29) 



n£Z\{0} n£N 

The (negative) Newtonian gravitational potential energy is given by 



U 



1 



■ h tt M . 



(4.30) 
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We may motivate this expression from two perspectives: geodesic motion or the red-shift 
of energy. 

From the geodesic motion point of view U is obtained by integration over the New- 
tonian force that is needed to bring the black hole from infinity. The Newtonian force is 
determined from the Geodesic equation 

x a + T« v x fl x u = 0, (4.31) 

that in the Newtonian limit (weak field limit) gives 

where hyj is given in fl3.4p . Comparing with Newton's second law we write the Newtonian 
gravitational field as 

E = ^Vh tt , (4.33) 

thereby motivating fl4.30|) . 

Alternatively, the energy red-shift is given by 



U = M - M « d t - 8 T ~ (l - V 3 ?* J d t 

~ (l - M = (l - y/1 + hu) M ~ -~ ht M . (4.34) 

Now we may substitute 

1 p d - 3 M 
n ~ 2(nL) d ^ 

and find 

r L - L d 4 - -L * V , j 1*?° , - - 3 > «" - 3 ) ^M„. (4.35) 

n G N 

Finally, substituting the expression for the rest mass ( [4.28 ) we recover the tension formula 

dm. 



5. Implications for the phase diagram 

In this section we shall translate the PN thermodynamic constants which we obtained into 
the phase diagram. The tension formula will dictate the slope of the small black hole 
branch. By extrapolation we shall attempt to learn about the phase transition region. 
Following the prediction for a critical dimension in this system we shall explore the 
dimensional dependence of the extrapolation and indeed we shall find an indication for 
Sorkin's critical dimension 

8 This extrapolation could have been done ever since the results were first found in 



13|.* 
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For the phase diagram it is convenient to define a dimensionless mass to serve as a 
control parameter 

G d M 

From the expression for the tension ( |4.25| ) we obtain an estimated Jl as a function of the 
dimensionless tension n for small black holes 

(d - 2)Q d . 2 

^ = 8vr(d-3) C (rf-3) (5 ' 1} 
where the dimensionless tension n is defined by 

U : ~ ~M' 



This quantity serves us as a good dimensionless order parameter for the phase diagram [10 
20], and on the uniform black-string branch it attains the constant value 



n„ = ^3. (5.2) 

In figure |] we give the mass \i against the tension n for d = 6 for both the analytic 
linear approximation ( |5.lD and the numerical results for the whole phase diagram (the black 
hole and the non-uniform black string) taken from |Q ] (where the results of [|TT] , |22|| are 
incorporated 9 ) . For the small black hole branch (lower right) we see excellent agreement. 

Next we would like to extrapolate and define the extrapolated intersection point of the 
linear approximation with the black hole branch 

»x := »BH{n st ) = 87r(d _3 )2c(d _ 3) , (5-3) 

(see also figure^]). Now we can compare fix with the Gregory-Laflamme (GL) critical mass 
fi cr it as a function of the dimension. In figure |5| we plot fix/^crit for various d. The critical 
mass was computed from Ugl-i the critical wave-number through 



(d - 3) ft d _3 



^ rit = - 16 7"° iff) • (5 - 4) 



d-4 



and values of kcL for various dimensions were obtained in |13|] and can be found in |2J 

For large d, kcL behaves as kci ~ Vd and the ratio fix/ ficrit tends to zero strongly. 
A similar observation appears in Q. We interpret that as an indication for a second order 
phase transition, where the black hole phase does not reach the Gregory-Laflamme region, 
but rather turns first into a stable non-uniform string phase, and only the latter joins into 
the GL point. Since |24j it is known that at d = 5 the transition is first order. Thus we 
find an indication for a critical dimension, which separates first and second order behavior 



and was indeed discovered by Sorkin to lie at D* = "13.5" |13j]. Interestingly, our graph 
shows a maximum in the range 12 < d < 14 which is the location of the critical dimension 
D*. 



9 We are grateful to T. Wiseman and H. Kudoh for sending us their data. 
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^ n/n st 



Figure 4: The analytic linear approximation line (red) compared to the numerical data points 
(green) of |Q for the complete phase diagram in d = 6 consisting of two branches: the black hole 
and the non- uniform black string. The vertical axis is proportional to the mass, and the horizontal 
axis to the tension (see text for exact definitions). For small black holes (lower right) we see 
excellent agreement. By extrapolation, the thick (blue) point defines /j,x, the intersection point of 
the extrapolated linear approximation with the line of uniform black strings, n = n sf . 



flfm ( |1=n st) 



l , , , , f d 

10 15 20 25 30 

Figure 5: The ratio Hx/^crit of extrapolated vs. actual Gregory-Laflamme points for various 
dimensions. For large d the ratio tends to zero strongly. We interpret that as an indication for 
a higher (second) order phase transition, and thus as an indication for Sorkin's critical dimension 
D* = "13.5". Interestingly, our graph shows a maximum around 12 < d < 14 which coincides with 
D*. 
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A. The Calculation of Pf(/ ($ z ) 2 dF d _i) 



We start by replacing $ by a sum over images 



Pf / (QvYdVd^ 



(A.l) 



Pf|(d-3) 2 ^ (d - 3) £ 

(m,n)G Z : 



dz 



„d-3 



dr- 



(z + mL) (z + nL) 



(r 2 + (z + m L) 2 ) 2 (r 2 + ( z + nL) 2 ) 2 



The summation here goes over all the pairs of mirror images of the black hole, and we 
transform it as follows 



. £ 

(m,n)G I 



dz 



where according to ([4.4|j4 



k=m— n£ I 



dz 



|fc=0 
00 

2 E 

k=l 



|fc=0 



fc=i 



+00 



(A.2) 



(A.3) 
(A.4) 



This transformation, replacing the integration over —k < z < k by integration over the 
covering space — 00 < z < 00, is justified since the integral is invariant under translations in 
the direction of the periodic coordinate z, and thus we can replace summation over images 
with fixed k = m — n by a summation of integrals over stripes (m — 1/2)L < z < (m+l/2)L 
which altogether exhaust the whole range of z in the covering space: for fixed k we have 
Yjm f~L/2 dz = J2m I(m-i/2)L dz ^,0, where l m .n is the integrand in (|Al|). Then we 
wrote the sum over k as a sum over natural numbers k € N (and multiplied by two) when 
k L is the distance of the black hole from its mirror image, and we separated the singular 
term which corresponds to k = 0. 

Now, we have to regularize just the singular term that comes from the self-interaction 

Pf (| dV d -^j = Pf (| {$>, Z 2 ) SI dVd-i) + j {®/) REG dV d ^. (A.5) 

This is done according to the idea of Hadamard's partie finie regularization introduced 
in section ||. Transforming the singular term to polar coordinates (p, x) an d imposing a 
cut-off at p = e we see that it is proportional to 



dp 

n d-2 



1 



oc 



T d-3 



namely 



Pf (/ Ms/ = 0. 



(A.6) 
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We turn to the regular part 

poo f'OO 

= 2p 2 {d - 3) (d-3) 2 y2 dz r d - 3 di 



(z + kL) 



feeN" 



(r 2 + z 2 ) 2 ( r 2 + (2 + A;.L) 2 )" 



— -(A-7) 



Integration by parts gives 



2 pi ^ (d - 3) £ f°d * / V 3 d r L__ . * f - 

£^/-co 7o (r 2 + z 2 )V ^ 



This integral can be rewritten in the form 



2 2(d-3) V- u_ 



d 2 



feel 



<% 2 



oo />oo 



b=kLJ-oo JO 



dz I r d 3 dr 



iZ + kL) B ■ (A, 
r 2 + (z + A:L) 2 ) 2 / 



1 1 

5^3 " ' d^F — ' (A.9) 

(r 2 + z 2 ) 2 (r 2 + (z + b) 2 ) 2 



transforming to polar coordinates and introducing a dimensionless variable u := £ leads us 
to 



= -2^-3)^ £ 



fceN 



a2_ 

9& 2 



1 



=feL 



lidii / d X ^ o;- (A.10) 

(u 2 + 2ncos(x) + 1) 2 



The integral over the angle x gives 

sin d ~ 3 (x) 



dx- 







(u 2 + 2ncos(x) + 1) 2 " 



1 



< U < 1 

— 1 < u < 00 



(A.ll) 



Finally, performing the 6-derivatives and integrating over u gives us the final result 



pfYy {®, z fdv d - 



2(d-3) 

.^-3-0^2^-4)^-3)^-3). 



(A.12) 



B. The next to leading correction to the "Archimedes" effect 

In ||] we computed the leading order for the inter-polar distance defined to be the proper 
distance between the "poles" of the black hole measured around the compact dimension 



L/2 

j poies = z / dz ^Jg zz , 
ha 



(B.l) 



where zjj denotes the location of the horizon. It is convenient to define the dimensionless 
quantities 



y: = 1 



^poles 

L 



PH_ 

L 



(B.2) 
(B.3) 
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where pn is the location of the horizon in isotropic coordinates, related to the Schwarzschild ra- 
dius, po through 

p Q d - 3 = Ap^ 3 . (B.4) 

In H we found 

y = 2I dV + o(r l ) , (B.5) 
where the constants Id are defined by 

^-i-/ 1 ((i + - d - 3 )^-i) 5 = 4lAv;F F(piy' k = d ~^ (R6) 

and are a monotonic function of d: I§ = 0, Iq = 0.6845, = 1. 

Here we shall compute further corrections for y, getting up to order d — 2 in r/. Let 
us start by identifying the necessary orders at each patch. In order to compute y at order 
r] s it is necessary to know the asymptotic metric up to order s, but in the near zone it is 
sufficient to know the metric up to order s — 1, due to the division by L in the definition 



(B.2). Thus for the leading result flB.5| ) it is sufficient to use the asymptotic zone metric up 
to order 1, which is the same as order 0, namely the flat metric, and in the near zone up to 
order 0, namely Schwarzschild. Here we will use our information about the metric in both 
zones up to order d — 3, which includes the Newtonian approximation in the asymptotic 
zone, and the first (monopole) correction in the near zone. Since the next order in the 
asymptotic zone always vanishes, our result will hold up to r] d ~ 2 . 10 

The way to compute L po i es is to pick some mid-point Z, divide the integration between 



the two zones 



2 



lpoies = / dz\j g y zz ' + I dz\ gi z J , (B.7) 



and confirm that the result is independent of the choice of mid-point. 
In the asymptotic zone we have (pj 



if ymp) = l + ^ + o(p^- 3) ) • (B.« 



Thus the contribution to L po \ cs / (2L) ( |B. 7| ) is 

1 fL/2 dz\Jgt ymp) 




L 



z 



„d-3 \ 

PO \ ,/r, ( 2(d-3) 



2(d-3) \z + nL\ d ~ 3 



O [Po 



10 In principle, we have enough information to get up to order 2(d — 3) — 1, since there are no corrections 
in the asymptotic zone up to that order and in the near zone there are known corrections from matching 



with the Newtonian potential (see subsection 3.3). However, we did not perform this computation since 
in 5d, the only order where we do not know yet the leading behavior since 1$ — 0, order 2(d — 3) — 1 = 3 
coincides with order d — 2 which we compute here, and thus such a higher order computation would add 
information only at higher dimensions where it is less needed. 
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+ 



Pt 3 



2 L 2(d-3)(d-A)L d - 3 



+o[p 2 (d - 3 \ 



oo 



^ (n + z) 

n=l v ' 



E 



Z/L 



+ 



J 1/2 



(B. 



where in the last equality we defined z := z/L. For 5 = 1/2 the expression in brackets in 
the last line vanishes. At the other boundary, z = Z/L, we should expand up to order z 
(since we need terms up to order L~( rf-2 ) in the overlap region. Note that here is quite easy 
to extend the result to order 2{d — 3) — 1, which we discussed above, simply by expanding 
further). Altogether we find 



L 



poles 

2L 



1 Z 



Po 



asymp 



2 L 2(d - 3)(d - 4) L d " 3 



d-4 



2(d-4)C(d-3)| 



. (B.10) 



We now turn to the near zone. Here at order (d — 3) we consider only the monopole 
correction, which leaves us with the Schwarzschild metric, with adjusted parameters. For 
matching purposes it is convenient to write the Schwarzschild metric in isotropic coordi- 
nates 



ds z 



1 - V 



dt 2 + (i + ^ (d P 2 c + p 2 c dn 2 d _ 2 ), 



where 



Ph\ 

Pc J 



d-3 



(B.ll) 



(B.12) 



To incorporate the correction we could have used ( |3.18| ) in a straightforward manner, 
but here we take an alternative route. In order to match we are free to adjust pn, rescale 
t and re-parametrize p c . Since we wish to retain the isotropic form we consider only a 
rescaling of p c . Matching of g zz is achieved by the transformation 



p c -> (1 + 5) p c , 

2(d-3) d-3 \LJ 



5 :-- 



(B.13) 



where $o is the Newtonian potential at the origin due to the images. Matching of gu 
requires also 

i-> (1 + (d-3)5)t . (B.14) 

One is free to change pn as well - this is a change of scheme, or reparametrization of the 
branch of solutions. Here it is convenient to choose pu to remain unchanged 



Sph = 



(B.15) 



(compare this scheme with the scheme of ( |3.18| ) which guarantees zero matching at PN 
order and in the current context would be expressed as pu — > Ph (1 + &)•) 
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We may now calculate the contribution of the near zone to L po \ es /(2L) 
-t . x 1 rZ 



d-3\ 2 /( rf - 3 ) 




(1 + *) 



dz 
Ph 



J PH 

p H /(Z(l+5)) 

f(l + t d-3)2/(d-3) 



(B.16) 



where in passing from the second to the third lines we used the adjustments ( B,13| ), in- 
cluding the adjustment in the location of the horizon induced by the rescaling of p c ; in the 
next line we performed a change of variables (1 + 8) z — > z and in passing to the fifth we 
changed to t := ph/z. Finally expanding the integral up to order p'jf 3 we find the total 
contribution from the near zone to be 



J poles 

2L 



near 



PH_ 
L 



Z(l + S) 

PH 



(d-3)(d-4) \Z(1 + S) 



Ph 



d-4 



(B.17) 



Adding up ( B.10 ) and ( |B.17|) we find that indeed all the Z dependence disappears 
to the prescribed order once we recall the definition of 5 ( |B.13| ) and the relation between 



Ph,Po ( B.4 ), and we are left with the final answer 



L 



poles 



2L 



--rjld + o 



f) 



d-2 



or equivalently 



y = 2I dV + o[{f 



d-2 



(B.18) 



(B.19) 



Summarizing, our result is expressed in terms of y, rj, Id whose definitions are given 
in ( |B.2| ,B.3,B.6). It turns out that using the scheme choice ( |B,15 ) our higher order result 
( B.19| ) is precisely our low order result (|B.5| ) . In particular we find that in 5d y vanishes up 
to order d — 2 = 3 in r\ (this is scheme independent) which is consistent with 5d numerical 
simulations [25| where y appears to be C(?? 4 ) (however, the result in ]7| is different). 
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